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ABSTRACT 


Cirrus cloud radiative and physical characteristics are determined 
using a combination of ground-based, aircraft, and satellite measurements 
taken as part of the FIRE Cirrus Intensive Field Observations (IFO) during 
October and November 1986. Lidar backscatter data are used to define cloud 
base, center, and top heights and the corresponding temperatures. 

Coincident GOES 4-km visible (0.65 pm) and 8-km infrared window (11.5 ^m) 
radiances are analyzed to determine cloud emittances and reflectances . 
Infrared optical depth is computed from the emittance results. Visible 
optical depth is derived from reflectance using a theoretical ice crystal 
scattering model and an empirical bidirectional reflectance model. No 
clouds with visible optical depths greater than 5 or infrared optical depths 
less than 0.1 were used in the analysis. 

Average cloud thickness ranged from 0.5 km to 8 km for the 71 
scenes. Mean vertical beam emittances derived from cloud-center 
temperatures were 0.62 for all scenes compared to 0.33 for the Case Study 
(October 27-28) reflecting the thinner clouds observed for the latter 
scenes. Relationships between cloud emittance, extinction coefficients, and 
temperature for the Case Study are very similar to those derived from 
earlier surface-based studies. The thicker clouds seen during the other IFO 
days yielded different results. Emittances derived using cloud- top 
temperature were ratioed to those determined from cloud-center temperature. 

A nearly linear relationship between these ratios and cloud- center 
temperature holds promise for determining actual cloud- top temperatures and 
cloud thicknesses from visible and infrared radiance pairs. 

An average visible scattering efficiency of 2.1 was found for this 
data set. The results reveal a significant dependence of scattering 


efficiency on cloud temperature. Values of mean scattering efficiency as 
high as 2.8 suggest the presence of small ice particles at temperatures 
below 235 K. Large uncertainties in the optical parameters due to cloud 
reflectance anisotropy and shading were found by analyzing data for various 
solar zenith angles and for simultaneous AVHRR data. The results highlight 
the need for additional study of cirrus cloud scattering and remote sensing. 



1. Introduction 


Accurate quantification of cirrus cloud properties from satellite 
measurements is particularly important to the understanding of the role of 
cirrus in climate change. The nonblackness of cirrus at thermal infrared 
wavelengths renders the interpretation of satellite data taken over cirrus 
more difficult than measurements over most water clouds. The International 
Satellite Cloud Climatology Project (ISCCP; see Schiffer and Rossow, 1983) 
is making an ambitious effort to derive daytime cirrus coverage, altitudes, 
and optical depths over the globe during a 5 -year period. The ISCCP 
analysis algorithm (Rossow et al. , 1988) relies entirely on bispectral data 
taken at visible (VIS, - 0.65 jim) and infrared (IR, - 11.5 pm) wavelengths. 
Although VIS-IR bispectral techniques have been suggested as feasible 
methods for determining bulk cirrus properties (e.g., Shenk and Curran, 

1973; Reynolds and Vonder Haar, 1977), there has been very little 
application of these techniques to real data prior to the ISCCP. 

The basic premise for using the bispectral approach is that the VIS 
extinction coefficient is related to the IR absorption coefficient. This 
relationship implies that the cloud VIS reflectance may be used to infer the 
cloud's IR emittance. Having a value for the clear-sky IR radiance, it is 
possible to correct the observed cloudy radiance for cloud emittance 
resulting in an estimate of the radiance emanating from some specified level 
in the cloud. The equivalent blackbody temperature of this level, usually 
the cloud center, is then converted to cloud altitude by means of a vertical 
sounding. The critical relationship ultimately required for this approach 
is the dependence of IR emittance on VIS reflectance through the IR and VIS 
optical depths. Since clouds scatter radiation anisotropically , this 
relationship is also influenced by the viewing and illumination conditions. 
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The ISCCP cirrus analysis (Rossow et al. , 1988) utilizes a combination of 
theoretical and empirical models to determine the cloud visible optical 
depth from the observed reflectance, the cloud emittance from the visible 
optical depth, and finally, the cloud- top temperature from the cloud 
emittance and the observed infrared radiance. The theoretical cloud model 
is a radiative transfer scheme which simulates the scattering and absorption 
of visible radiation by water droplets with an effective radius of 10 p m • 

For water droplets of this size, the ratio of VIS extinction to infrared 
absorption is - 2.7. An analysis of coincident satellite and lidar data by 
Platt el al. (1980; hereafter, PRA) and theoretical calculations employing 
cylinders (Platt, 1979) suggest that this ratio is approximately equal to 
2.0 for cirrus. The ISCCP algorithm utilizes the latter value to provide a 
link between the water droplet model and actual cirrus clouds . 

Cirrus clouds are primarily composed of ice crystals with various shapes 
having maximum dimensions ranging from about 20 pm to 2000 pm (e.g., 

Heyms field and Platt, 1984). The scattering properties of hexagonal ice 
crystals differ considerably from spherical particles (Liou, 1986). Because 
of the complexities involved in computing scattering by hexagonal solids , 
cylindrical columns have been used to approximate hexagonal crystals in 
radiative transfer calculations (e.g., Liou, 1973). More recently, however, 
Takano and Liou (1989a) have solved the radiative transfer equations for 
randomly oriented hexagonal plates and columns. Their results are the most 
realistic to date in that they reproduce certain well-known cirrus optical 
phenomena . 

Absorption plays the dominant role in IR extinction in cirrus clouds. 

Some theoretical investigations (Liou and Wittman, 1979; Stephens, 1980), 
however, have shown that scattering effects may also be significant at IR 
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optical depths greater than - 0.1. An IR radiance measured by a satellite 
over cirrus clouds, therefore, is the product of both absorption and 
scattering processes in the cloud as well as the transmission of radiation 
from below the cloud (Platt and Stephens, 1980). These effects may 
complicate the estimation of emittance from the observations. The emittance 
derived from the satellite data is an effective beam emittance. Some 
corrections for scattering effects may be needed to obtain the absorption 
beam emittance which determines the cloud IR optical depth. 

Empirical studies have also shed some light on the VIS reflectance - IR 
emittance relationship. Platt (1973) developed techniques for deriving 
cloud visible and infrared properties from a ground-based lidar and an 
upward- looking infrared radiometer. The backscattered intensities measured 
with the lidar are used to define cloud base and top heights. Cloud 
emittance was derived from the observed downwelling IR radiance . Platt and 
Dilley (1979) presented emittance results from a set of observations taken 
over Australia. PRA used lidar and satellite VIS-IR data to estimate the 
dependence of cloud beam emittance on cloud VIS reflectance for a limited 
set of viewing and illumination conditions over Colorado. Their results are 
more consistent with theoretical scattering from ice cylinders than with 
scattering from ice spheres. Aircraft radiometric measurements taken over 
New Mexico (Paltridge and Platt, 1981) have also been used to determine the 
radiative characteristics of cirrus clouds as related to the cloud ice water 
path. Those results provide further evidence that real clouds scatter more 
like cylinders than spheres. Platt (1983) combined the results from 
previous studies and used them to explain the characteristics of two- 
dimensional bispectral histograms of VIS-IR data observed from a 
geostationary satellite. Theoretical calculations of reflectance and 
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emittance for typical cirrus clouds were consistent with the satellite data 
taken over areas of suspected cirrus clouds . While that study provided 
encouragement for using a bispectral approach to retrieving cirrus 
properties from bispectral data, it also highlighted some of the 
difficulties which are likely to be encountered with such a technique. 

Platt and Dilley (1984) used lidar and solar radiation measurements to 
measure part of the single-scattering phase function of real cirrus clouds. 
Their results fell within the range of laboratory measurements and 
theoretical calculations for hexagonal crystals. An analysis of a large 
sample of ground-based lidar and infrared data taken over Australia (Platt 
et al., 1987) showed that the average emissivity of cirrus clouds is 
primarily a function of the midcloud temperature. Though fraught with 
significant uncertainties, that study also indicated that the theoretical 
value of the ratio of visible extinction to infrared absorption for cirrus 
clouds may be too low. 

From these previous studies , it appears that 

(1) cirrus cloud scattering properties are similar to those of 
hexagonal crystals resulting in reflectance patterns which are 
unlike those from spheres; 

(2) scattering of IR radiation may be important in determinations of 
IR optical depths ; and 

(3) the ratio of VIS extinction to IR absorption coefficients is 
between - 1.8 and 4.0. 

The full impact of these results on using a VIS-IR bispectral method for 
retrieving cirrus properties is unknown. Differences between ice crystal 
and water droplet bidirectional reflectance patterns will introduce errors 
into the retrieved VIS optical depth. Use of an effective beam emittance 
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with a theoretical model which assumes absorption only may affect the 
emittance estimation. Finally, uncertainties in the extinction ratio 
(scattering efficiency) may cause significant errors in the estimation of IR 
optical depth. 

In this paper, the relationship between VIS reflectance and IR emittance 
is examined using data taken during the First ISCCP Regional Experiment 
(FIRE) Cirrus Intensive Field Observations (IFO, see Starr, 1987). Ground- 
based and aircraft lidars are used to define the vertical locations of the 
cirrus clouds while satellites provide measurements of VIS and IR radiances. 
Both VIS and IR optical depths are computed from the reflectance and 
emittance data covering a range of solar zenith angles missed in previous 
studies. These relationships are derived to provide a means for the 
application of a bispectral cirrus parameter retrieval algorithm over the 
FIRE IFO region. Results are presented for the entire IFO period with 
emphasis on October 27 and 28, 1986. The data presented here also 
constitute an initial source for developing cirrus bidirectional reflectance 
models and may be used to help validate the models employed in the ISCCP 
algorithm. 
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2. DATA 

a. Lidar measurements 

Lidar backscatter data were taken from four different sources- -three 
surface and one airborne. The lidars and their operating systems and data 
products have been described elsewhere. Thus, only a brief description of 
sources and their uses in this study are given here. 

The University of Utah mobile polarization lidar (see Sassen et al. , 

1989) was located at Wausau, Wisconsin (WAU, 45.0°N, 89.7°W). The NASA 
Langley ground lidar (Alvarez et al . , 1989) was situated at Ft. McCoy, 
Wisconsin (FMC , 43.9°N, 90.8°W), while the University of Wisconsin high 
spectral resolution lidar (Grund and Eloranta, 1989) was in Madison (MAD, 
43.1°N, 89.4°W). These ground systems acquired nearly continuous lidar 
backscatter profiles during the cirrus days of the IFO with especially good 
coverage during the case study period. The lidar returns are used to define 
the cloud base and physical thickness. Under conditions of small attenu- 
ation and constant backscatter phase function, the backscatter intensity 
profiles indicate the vertical distribution of cloud extinction. A time 
series of these lidar returns shown in Fig. 1 define the outlines of the 
cirrus clouds as they passed over FMC during the afternoon of October 28 . 
Solid black areas define the most intense cloud backscatter. Grey denotes 
less backscatter and white indicates no cloud. Vertical stripes represent 
missing data. Cloud-top altitude is fairly constant at - 11 km. Cloud base 
changes from - 8 km to 10 km approximately every half hour. In the morning, 
cloud base was observed at — 7 km, while cloud top varied between 8 km and 
11 km. A similar variation is also seen in the WAU estimated volume 
backscatter coefficients, fi*, shown in Fig. 2 for the morning of October 28. 
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Darker portions of the plot correspond to higher values of 0 * • The 

particle backscattering efficiency depends on cloud particle shape and 
phase. Further details of the lidar returns are reported in the cited 
references . 

Three parameters are derived from plots like those in Figs . 1 and 2 by 
averaging the data within + 15 min of the GMT (Greenwich Mean Time) half 
hour plus 5 minutes. All times, however, will be given here to the nearest 
half hour. Cloud- top altitude, z , and cloud-base altitude, z b> are defined 

as the average altitudes of the highest and lowest nonclear-air backscatter 
returns , respectively. Similarly, the cloud thickness is h ■ z^ - . 

Mean cloud height (approximately cloud center height), z c> is the 

backscatter- intensity weighted average height of the cloud. It corresponds 

roughly to the altitude below which 50% of the lidar backscatter is 

accumulated. These parameters were estimated graphically for the FMC and 

WAU sites, while a computer analysis was applied to the MAD results. The 

value of z for MAD corresponds to the midpoint in optical thickness 
c 

independent of cloud attenuation (see Grund and Eloranta, 1989). Since the 
clouds are advecting over the fixed surface sites, the averaged lidar data 
correspond to a thin vertical cross-section taken out of some cloud volume. 
It is assumed that the cross -section- averaged data represent the mean 
conditions of that volume. 

These same parameters were also derived from the down- looking lidar 
backscatter plots reported by Spinhime et al. (1988) for selected flight 
tracks of the high-flying, NASA ER-2 aircraft over the IFO area. Shorter 
time averages were used since the plane's motion greatly increased the 
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cirrus advection rates relative to the lidar. In some instances, the clouds 
were too thick for complete penetration by the ground-based lidars. To 
determine these occurrences, the cloud altitudes estimated from the ground 
were compared to those determined from the nearest aircraft flight. On most 
days, there was good agreement between the surface and airborne lidars. The 
thick clouds observed on October 22 required use of the aircraft lidar to 
obtain a good estimate of z^. At other times when no aircraft data were 

available for comparison, a different approach was used to estimate z t (see 

Sec. 3a). 

b. Temperature data 

Soundings from Green Bay, Wisconsin, were used to determine temperature- 
height relationships for all of the data. Linear interpolation was used to 
estimate half-hourly soundings from the 6 -hourly data . Cloud- top 
temperature, T , corresponds to z on the soundings. Mean cloud 

temperature, T c> is found from z q . Surface temperatures taken every 6 

hours at MAD, WAU, and Lone Rock, Wisconsin (Hahn et al. t 1988), and 
occasionally at the FMC site were used to supplement the clear- sky 
temperatures derived from the satellite data (see Sec. 3b). 

c. Satellite radiances 

Half-hourly VIS and IR data from the Geostationary Operational Satellite 
(GOES) were resolved to 4-km pixels by averaging and replication of 1-km VIS 
and 4-km x 8-km IR pixels, respectively, for most times. When the full- 
resolution data were not available, 8-km pixels derived by pixel averaging 
(VIS) and sampling (IR) were used. The pixels were navigated to + 8 km as 
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in Minnis and Harrison (1984a). Pixel pairs were extracted for areas within 
0.35° longitude and 0.25° latitude of the surface lidar sites. These areas 

correspond to cloud advection at wind speeds of — 30 ms for a half-hour 

interval. Strips of pixels, 4 (2) pixels wide, corresponding to the wind 

vector at z centered over the surface site were then taken from the 4 -km 
c 

(8-km) resolution areas. The satellite data corresponding to the ER-2 
flight track were extracted for a strip centered on the flight track. This 
strip is approximately 30 km wide and 100 km long. 

Two-dimensional histograms of the VIS and IR radiance pixel pairs were 
formed from the strips of extracted data. The VIS data, stored as counts, D 
(where 0 < D < 63), in the histogram, were converted to radiance and VIS 
reflectance, p, using the calibration of Whitlock et al. (1989). The raw IR 
data are given as equivalent blackbody temperatures, T, and are converted to 
radiance with the Planck function, B(T), evaluated at 11.5 p.m. Note that 
all radiometric quantities discussed in this study are spectral quantities , 
either VIS or IR, so no subscripts denoting spectral dependence are given. 

A similar set of histograms was derived from the NOAA-9 Advanced Very 
High Resolution Radiometer (AVHRR) Global Area Coverage (GAC) 4-km data. 

The NOAA-9 daytime pass over the areas occurred between 2000 and 2130 GMT 
(- 1430 local time) . Radiances from channels 4 and 5 were averaged to 
produce an IR radiance close to that of the GOES. The AVHRR VIS data were 
taken from channel 1 and converted to reflectance using the calibrations of 
Whitlock et al. (1989, personal communication). 

Solar zenith, satellite zenith, and relative azimuth angles, 8 0 , 6, and 

respectively, were computed for each set of measurements. The only GOES 
available at that time was GOES -6 located over the Equator at approximately 
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97°W. This location yields a value of 6 - 52° for the ground sites. 
Relative azimuth angles were confined to the backscattering hemisphere for 
GOES. The orbital precession of the NOAA-9 produces a variation in 6 
between 0° and 70° over a given site every 4-5 days. Because of its cross - 
track scan pattern, the AVHRR views a given site from relative azimuth 
angles which lie in both the forward and backward hemisphere at a nearly 
constant skew to the solar plane. The skew depends on both the solar and 
site latitudes. Consequently, a range of solar zenith angles with 
relatively constant viewing angles is covered using the GOES, while a range 
of viewing angles at a constant 8 0 is covered with the AVHRR. 

Days when cirrus were observed over the IFO without substantial low-cloud 
interference are referred to as cirrus days. The times used in this study 
when lidar and satellite data coincided with cirrus are listed in Table 1. 


3. Methodology 

a. Emittance 

The effective cloud beam emittance is given here as 

e b (f) - [B(T) - B(T g )]/[B(T z ) - B(T s >], (1) 

where T is the temperature at some altitude z corresponding to the 

cloud. The mean clear- sky equivalent blackbody temperature over the area of 
interest, T , has a weak dependence on 6 . The effective beam emittance 

Includes the effects of scattering as noted earlier. 

Cloud beam emittance is calculated twice for each set of lidar-radiance 

data using T - T and T - T . The former value, which corresponds to 
z c z t 
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the quantity used in most previous studies (e.g., PRA) , may be more 

representative of the actual radiating part of the cloud. It is unlikely 

that a value of c. - 1 will be measured at a useful satellite zenith angle 
b 

using T because of the low density of particles in the upper portions of 

the cloud. On the other hand, T c may be too great for some thick cirrus 

clouds. In those cases, > 1. Although effective beam emittances greater 

than unity may be possible due to scattering enhancements of the upward 
radiance (Platt and Stephens, 1980), the uncertainties in z q for thick 

clouds preclude any definitive measurements of > 1. Thus, if > 1, 

T is decreased until the average value of e, for a given reflectance is 
c D 

less than or equal to one. 

During initial processing of the data, it was determined that the maximum 
emittance found using T was - 0.86, except for those cases with cloud 

cover too dense for complete penetration of the lidar beam. To identify and 
correct the exceptions, a new estimate of T t was computed whenever T c 

was adjusted as explained above. This new estimate, 

T' - + [ B (T) - B(T ) ] / 0.86]}, 

t S 

was then compared to T . If T£ < - 3K. then T fc - T^. The result was 

then compared to the tropopause temperature. If T t was colder, it was set 

equal to the tropopause temperature. The value of was then adjusted to 

correspond to the value of T. . 

It is also assumed that 

< b - 1 - exp(-T e //i) , 


( 2 ) 
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where r is the effective IR optical depth and n — cos0 . Based on the 
e 

results of Platt and Stephens (1980) , it is expected that the viewing zenith 
Angle dependence of will not depart significantly from (2). Values of 

effective beam emittance derived with AVHRR data may be adjusted to the GOES 

viewing zenith angle with this relationship. 

The vertical effective beam emittance from (2) is 

c - 1 - exp ( - r ) . (3) 

a ° 

It is assumed here that scattering effects are negligible in the upwelling 
direction. Thus, r & is equivalent to the IR absorption optical depth and 

is equal to the vertical emittance. This assumption is consistent with 

the inability to confidently obtain £ b > 1 with the current data set. 

b . Clear -skv temperature and reflectance 

The clear-sky temperature is estimated in several different ways. The 
first order estimate is taken from the initial results of Minnis et al. 
(1989) which applies the techniques of Minnis et al. (1987) to 0.5° regions 
within the greater IFO area (42°N - 47°N, 87°W - 92°W). That approach sets 
a VIS threshold - 2 counts above the clear-sky reflectance, or count, 

D . All pixels which are darker than this threshold and have T > T max - 3 

K are considered to be clear. Their average is T cg . The surface air 

temperatures, T , are also taken from nearby ground stations. A rough 
S 

correction is applied to these temperatures to adjust for atmospheric 
attenuation and the difference between the temperature of the surface skin 
and the air at shelter height. The resulting estimate of clear-sky 
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temperature is T . An example of the relationship between T cg and 

is shown in Fig. 3 for a region including EMC during October 27, 1986. 

Note that the difference between T and T constitutes the correction to 

c s § 

determine T from T . It varies with local time from positive values in 
a g 

the morning to negative values in the afternoon. This variation is typical 

of land surfaces (e.g., Minnis and Harrison, 1984a). 

The initial value of T for the ground site is compared to the values 

cs 

for the surrounding regions to check for cirrus contamination of the clear- 
sky temperature. If T cg for the ground site is not within 2 K of the 

warmest value found for nearby regions , then its value is reset to that of 
the warmest value. Normally, T g - T cs - If extensive overcast prevents the 

calculation of T , T - T . 

css a 

The AVHRR clear-sky temperatures are first estimated using T cg derived 

from the AVHRR data. If that result is lower than the corresponding GOES 

estimate of T by 2K, then the AVHRR value of T is reset to the GOES 
s cs 

value. Since the value of T g should vary by only 1 or 2 K over the range 

of 8 used here, limb -darkening effects are assumed to be negligible in the 
AVHRR analysis. 

Values of clear-sky reflectance, p g , were computed for each region using 

the 0.01° clear-sky albedo, a g , map of the IFO area constructed by Minnis et 

\ 

al. (1989) from GOES data at each half hour. Clear-sky reflectance over any 
latitude, A, and longitude, 4 >, of the grid at time, t,is estimated as 
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where x g is the anisotropic reflectance factor with values given by the 
model of Minnis and Harrison (1984b). The value of 0 O varies by a few 
degrees over the IFO time period, while the values of a g were normalized 
to a single value of solar zenith angle designated To account for 

these variations, c* s (t,0 o ) — a g (t , 8 0 )Mo/Vo t > where /* 0 — cos0 o and /*o t - 
cos0 Ot - The clear-sky diffuse albedo is 

at sd - J a s (8 o)A»o d Mo / J Po d Mo - 

integrated over - 0,1. The value of a gd is estimated to be equal to 

the value of c* s (57 0 ) in this study. 



The transmittance of the air above the cloud is 

T - exp[-k u(1/m 0 + 1/70], 

& o 

where k - 0.085 - 0.00052u (Rossow et al., 1988) and u is the ozone 
o 

abundance in cm-STP. The value used here, u - 0.32 cm-STP, is the average 
of the midlatitude winter and summer standard atmospheres above 10 km from 
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McClatchey et al. (1973). PRA implicitly assumed that - 1, although 

there is significant ozone absorption in the Chappius bands. 

The transmittance of the cloud to direct solar radiation at 9 0 is 

T c - exp ( ' r v / 2/i 0 ) , (6) 

(see PRA). Similarly, the direct transmittance from the surface through the 
cloud along the satellite line of sight is 
T u - exp(-r v / 2/i) . 

The visible optical depth is reduced by a factor of two for the direct 
transmittance because at least half of the radiation scattered out of the 
beam is actually diffracted in the forward direction (Takano and Liou, 

1989a). Clear-sky reflectance along the satellite line of sight is p g , and 

a is the effective clear-sky albedo to diffuse radiation directly below 
sd 

the cloud. Due to the relative homogeneity of clear-sky reflectance over 
the IFO region, it is assumed that a gd and p g may be computed from the 

same data. The albedo of the cloud to diffuse radiation is q^. 

This model assumes that all of the ozone absorption occurs above the 
cloud (first term) and that all Rayliegh and aerosol scattering is confined 
to the layers below the cloud. The second term in (5) accounts for direct 
solar radiation which passes through the cloud, reflects from the surface, 
and passes back through the cloud in the direction of the satellite. The 
third term accounts for the radiation which passes down through the cloud 
via multiple scattering, reflects from the surface (which now acts as a 
Lambertian surface) below the cloud, and returns through the cloud scattered 
in the direction of the satellite. 
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In addition to values for the clear -sky terms, solution of (5) for p c 

requires specification of r and x c - Here, the value of the VIS optical 

depth is estimated by iteration on (5) using a log-linear interpolation of 
the relationships between p 0 and c« c for randomly oriented hexagonal 

columns (length, 125 pm; width, 50 pm) in Fig. 4 of Takano and Liou (1989b). 
Similar interpolations are used to estimate <* d ( r v ) . where 

«d<v> - J 0 °d ( V Mo)dMo / lo 

Given a measurement, (5) may be solved for an initial guess of cloud albedo 
by assuming that T c “ T u “ 1 ‘ Q c - A value for T v is determined from 

this initial guess using the theoretical data. It is then used in (5) to 

solve for a . The initial guess and resulting value of a are compared, 
c c 

If the absolute difference is greater than 0.001, the iteration proceeds 

using the new value of as the guess. A limit of 20 iterations is 

imposed, although fewer than 5 iterations are required to achieve a 
difference of ± 0.001 between the guess and the computed value of c* c for 

most cases. Since <* c may be greater than zero, even for p < p g , t* c is 

set to 0.001 for initial guesses which are less than or equal to zero. If 

r < 0 a < 0.001, or a < 0.001 after any iteration, it is assumed that 
v ’ d c 

r is indeterminate and the data are not used, 
v 

A value for x c is needed to determine c« c from p c » x c depends on 

r v and the cloud microphysics. No models of x c are currently available 
for ice clouds in terms of r . Because of favorable angles PRA were able 
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to assume thflt 1 . Most empirical and theoretical bidirectional 

reflectance models for cloudy scenes (e.g., Suttles et al., 1988), however, 
reveal a systematic decrease in x c with 6^ for the angles used in this 

study. The cloudy scene bidirectional reflectance model developed by Minnis 
and Harrison (1984b) is used to estimate x c - That model was derived for 

top-of-the-atmosphere (TOA) broadband shortwave reflectance from clouds of 
all types. It is assumed that the broadband reflectance patterns are the 
same as those in the VIS regime. This model's reflectance anisotropy is 
similar to other empirical and theoretical models (Stuhlmann et al . , 1985). 
The inclusion of all cloud types in its derivation should produce a 
reflectance pattern between those for ice and liquid water clouds. 

Examples of the variations of p as a function of a c computed with (5) 

are shown in Fig. 4a for WAU at 1400 and 1800 GMT. The corresponding values 
of p s are 0.13 and 0.16, respectively, with a gd - 0.11. The relationship 

between cloud albedo and optical depth depends on the solar zenith angle 
resulting in a divergence of the curves at larger values of a Additional 

calculations were performed for p 0 - 0.4, x c “ 1. and P s “ 0.04, 0.08, 

0.16, 0.32, and 0.64 with a gd - 0.92p g . The results are shown in Fig. 4b. 
For relatively dark surfaces, the measured reflectance increases mono- 
tonically with a c< Over bright surfaces, the reflectance actually 

decreases first, then increases at greater cloud albedos. A dark surface 
contributes little to the upward flux, while the opposite is true for bright 
surfaces . The increased reflectance of a cloud to diffuse radiation tends 
of the reflected radiation from the surface. For thin clouds 


to trap some 



18 


over bright surfaces, this effect can result in a reflectance which is lower 
than that observed for the clear-sky case. Additional calculations for 
other solar zenith angles indicate that the insertion of the thin cirrus 
cloud above the bright surface decreases the albedo at lower sun angles 
making the discrimination of clear and cloudy skies more difficult over 
deserts and other bright scenes. It is possible that this effect may also 
be important at certain angles and values of x c over dark surfaces . 

d. Scatterine-to-absorption ratio 

The mean VIS scattering to IR absorption optical depth ratio 
N 

r - 2 (r / r ) / N (7) 

i-1 i 

is computed for each value of p . This parameter is equivalent to the 
scattering efficiency factor. Only one temperature is used to compute 

for a given data set since only one average cloud height is derived for each 
time. Changes in the actual cloud height and thickness within the scene 
(e.g., Fig. 1) tend to introduce variations in for a given reflectance. 

A mean value of r is computed for each cloud reflectance value to minimize 

the effects of cloud height variability. 

If it is assumed that a c a for small optical depths, then 

« b - 1 - exp(-ka c //i) , (8) 

where k is the constant of proportionality. This relationship is 
primarily used here to filter low clouds. Values for k are determined in 
the following manner. For a given site and time, « b is computed with (1) 

for each pixel with D > D and T < T - 

s s 


3 K. The mean emittance is then 
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computed for each value of D corresponding to a value of c* c . Equation 

(8) is solved for k using each albedo -emittance pair. An average value of 
k is then computed for a specified subset of the data. This procedure is 
applied using both T z - and T z - T^. 

e. Low and thick cloud filtering 

The presence of low clouds contaminates the data since the analysis is 
predicated on the presence of only one cloud layer (it may contain two or 
more distinct sublayers). The cloud heights are defined for the upper-level 
cloud (generally, z fe > 4 km) detected by the lidar. In some instances, the 

occurrence of low clouds was either noted by the lidar operators or evident 

in the backscatter intensity plots. In other cases, the visual observations 

were not recorded or the low clouds were located within the satellite - 

defined area around the site, but beyond the observer's horizon. The data 

sets containing low clouds must either be eliminated or the pixels 

contaminated by low clouds must be removed. Elimination of all data sets 

containing low clouds would substantially reduce the number of samples. 

To filter the data, it is assumed that low clouds are generally brighter 

than cirrus, but their emittances depend on albedo, as in (8), in the same 

manner as cirrus. Thus, if a maximum cloud albedo is known for a given 

cirrus emittance, then it may be assumed that any pixels brighter than that 

maximum contain some low clouds. The values of those maximum albeods, o , 

nioX 

must be estimated here. To determine those maximum values, the histograms 
were first filtered manually to obtain a set containing no low clouds. The 
maximum albedo for each temperature (emittance) was extracted for every 
histogram at a given GMT. Data taken on October 22, 1986, were used in the 
analysis. All maxima for a given GMT were compared to determine the 
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greatest values for that time. Very little dependence on GMT was found in 
this comparison. Therefore, the maxima from all hours were combined. The 
resultant values were then used to solve (8) for k. The results were 
averaged to obtain the coefficient, k^, which is used to define the 

maximum cloud albedo for a given emittance. The results of this process 

yield k - 2.4 for T - T and k - 2.0 for T - T Thus, for a 
J max z c max ^ L 

given time, (8) may be solved for a max using k^ and the emittance. 

Anv pixel with a > a + 0.025 is assumed to contain low clouds and is 
J K c max 

eliminated from further processing. 

An example of this filtering is shown in Fig. 5 for the data over FMC at 
1930 GMT, October 28, when low clouds were noted in the observer's log. 
Crosses represent individual pixels. The solid line represents Q ma x 

computed using k The vertical dashed line, o - 0.20, is used at low 

r ° max max 

reflectances because of large uncertainties at those reflectance levels. 

The filtering process eliminates the pixels which clearly contain low clouds 

(i.e., those to the right of the line). Some pixels containing both low and 

high clouds, however, may remain in the data. 

All cases having > 5.0 were also removed from the data set. This 

arbitrary limit was imposed to minimize the uncertainties in r , r^, and 

subsequent values of r. The error in derived optical depth increases with 
increasing emittance. The value of T c is also more uncertain for thick 

clouds than thinner clouds. Thus, errors in emittance are greater for these 
thick clouds. Their removal does not significantly affect the trends in the 
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results. Inclusion of the thick clouds, however, greatly increases the 
noise in the data. No case study scenes were affected by this filtering. 

f . Partially cloud-filled pixels 

The use of (7) or (8) implicitly assumes that the pixel is completely 
filled with cloud. While it is likely that some of the pixels are partially 
clear, definitive means for determining which are overcast and which are 
partly cloudy are unavailable. Platt (1983) proposed a bispectral method to 
discriminate partly cloudy pixels from those containing thin clouds. His 
technique requires uniform blocks of constant height and thickness for 
detection of partly cloudy pixels. Detection of thin clouds allows variable 
thicknesses, but not heights. As seen in Figs. 1 and 2, real cirrus clouds 
are anything but uniform and constant in altitude and depth. An actual 
cloud field may also contain broken clouds with variable optical depths. 

The nonuniformity of such real clouds blurs the distinction between the 
linear relationship for partly cloudy pixels and the exponential dependency 
for variable thickness overcast pixels. For consistency with the ISCCP 
algorithm, it is assumed here that the pixels colder than T |nay - 3 K are 

overcast. The data are not used, however, if there is evidence of any 
significant breaks in the ceiling from the lidar displays. Very few breaks 
were seen over the sites on October 28. The cirrus observed on October 27 
were mostly scattered leading to the selection of only two scenes for that 
day. 
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4. Results 

a. GOES emittances for midcl oud temperature 

All results discussed in this section are based on T z - T c in (1). 

Examples of the two-dimensional histograms used in this analysis are shown 
in Figs. 6 a and 6 b for 1500 GMT over FMC . The latter represents a cirrus 
case (see Fig. 1) on October 28, while the former, taken during the previous 
day , ^5 typical of clear conditions. Maximum clear-sky reflectance for this 
hour is denoted with the dashed line in Fig. 6 a. Some of the cold, 
apparently cloudy pixels in Fig. 6 b are no brighter than the clear pixels in 
Fig. 6 a. Moreover, some of these pixels are actually darker than the cloud- 
free pixels. Depending on x c and r^, some of the cold, dim pixels yield 

a positive value of <» c in the solution of (5) . Those pixels with “ c < 0 

or T < T - 3K are, hereafter, referred to as "dark" pixels. They are not 
s 

used to solve either (7) or ( 8 ). Their impact and origins are discussed in 
a later section. The cloud emittances are plotted in Fig. 7a against the 
measured reflectances for the case in Fig. 6 b. Eliminating the "dark" 
pixels and applying (5) to the data in Fig. 7a yields the pixels values 
shown in Fig. 7b for cloud albedo. Averaging the emittances for a given 
albedo produces the mean values plotted in Fig. 7c. The solid line 
represents ( 8 ) using the average value of k - 5.1. The mean effective 

emittance and visible optical depths for the data used in Fig. 7c are « b - 

0.38 and - 0.59, respectively, yielding r - 1.73. 

The data from FMC in Fig. 7c for this hour are compared to those from MAD 
and WAU shown in Fig. 8 . Apparently, the clouds over WAU are much denser 
than those over MAD while the MAD observations are similar to those over 
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EMC. Values of T c differed by only 1 K among the sites, while T t ranged 

from 225 K at WAU to 217 K at MAD. Depolarization ratios derived from the 
lidar returns indicated intermittent liquid layers during the morning of 
October 28, especially at ~ 1500 GMT. The data from all three sites were 
combined, averaged, and fit with (8) yielding k — 5.6. Averaging the 
optical depth ratios for all three sites yields r - 2.01. The scatter in 
the means between the sites at a given hour is of the same order as that for 
different hours at the same site as seen in Fig. 9 for FMC at 3 times during 
October 28. 

A summary of the results for the Case Study (October 27-28) , hereafter 
referred to as CS, are given in Table 2. Cloud optical depths were much 
greater over WAU than over the other sites. The percentage of dark pixels 
was higher over FMC and MAD than over WAU. Cloud-top heights range from 9.5 
to 11.0 km at all three sites. The cloud-center temperatures vary by about 
25 K. 

Due to dropouts, the only data available for 1500 GMT during the primary 
IFO cirrus days occurred on October 28. Data from other days were available 
for most of the afternoon hours. For example, the cloud albedo -emittance 
pairs taken at 2000 GMT (Fig. 10) cover the full range of emittances . All 
of the data used in the analysis for the CS are shown in Fig. 11a. The 
clouds over the area during October 27 and 28 were generally much thinner 
with lower emittances than most of those observed during the other days . 
Combining data for all of the times given in Table 1 produces a large number 
of samples for « b < 0.5 and e b > 0.8. No data are found for « b < 0.08 

since no clouds are recognized if T > T g - 


3 K. 
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Table 3 summarizes the values of r v and r derived from all of the IFO 

data (Fig. lib) and from CS data only (Fig. 11a) for each relevant GMT. The 
average scattering angle, 8, between the sun, satellite, and scene is also 
listed in Table 3. Visible optical depths observed during the CS are less 
than half of those observed for all of the IFO cirrus days. Despite these 
and other differences, two similarities are quite evident. For both data 
sets, the scattering efficiency appears to increase with decreasing 0 Q and 

increasing 0. At high values of 0 O , r is well below the expected limit 

of 2. The average values of r are also very close, 2.07 and 2.08, for the 
CS and IFO, respectively. 

The various values of VIS and IR optical depths derived from Fig. 11 were 
averaged for discrete levels of r . Figure 12 shows the variation of mean 

r with IR optical depth. The standard deviations are denoted with the 
vertical lines. The CS data (squares) vary from 1.62 to 2.55, while the IFO 
data vary range from 1.60 to 2.31. There is good consistency between the 
IFO (circles) and CS results for < 1.0. Clouds with r ^ > 1.0 were 

undersampled during the CS (Fig. 11a). The trend of decreasing values of r 
with for < 1.0 may be due to dark pixel effects (see section 5c). 

Average scattering efficiencies were also computed for discrete levels of 
T . They are plotted in Fig. 13. The value of r is close to 2.0 for T c 

> 235K. Maximum values of r, approximately 2.8 for both data sets, were 

found for T « 225 K. The lower value of r for T - 217 K is derived 
c c 

from only five cases which were mostly taken during the late afternoon. 


25 


The mean vertical beam emittance is given as a function of temperature in 
Fig. 14. Although there is a general increase in with increasing cloud 

temperature, the values derived for the CS are generally lower than the IFO 
means. Cloud thickness versus T c is shown in Fig. 15. Maximum average 

thickness is found at T ~ 232 K for the IFO and at T s 241 K for the CS . 

c c 

The CS values of h increase nearly monotonically with T c - The IFO 

results show a tendency towards thinner clouds for T fi < 220 K and for T c > 

240 K. The thickest clouds were observed during October 22 and November 1. 
Cloud IR volume absorption coefficient, a - r / h, is plotted against T 

in Fig. 16. The results reveal a trend of increasing extinction with T c 

in the IFO data (correlation coefficient of 0.71). The average value of 

is 0.20 for the CS data. 

b. GOES emittances for cloud- top temperature 

The analyses discussed above were also performed for the GOES -derived 
emittances for T - T fc . A plot of all of the mean cloud emittance -albedo 

pairs is shown in Fig. 17. The largest concentrations of data are found for 

e, (T ) < 0.8. In general, o is greater for a given value of e, than it 
b t c D 

is in Fig. lib. Cloud- top emittance tends to plateau at « 0.89, while 

cloud- center emittance appears to level at * 0.98. The lower emittances 

lead to diminished values of r g and greater values of r relative to 

those derived for T . 

c 
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The emittance ratio, X ( - « b (T t ) / « b (T c ), was computed for discrete 

intervals of T . Mean values and standard deviations of these ratios are 
c 

shown in Fig. 18. The emittance ratio increases almost linearly with 
decreasing cloud center temperature. Standard deviations about a given mean 
ratio are less than 0.1. The emittance ratio is close to unity for < 

215 K. 


c. GOES/ER-2 cloud parameters 

Table 4 lists the data and derived parameter values for the five cirrus 
matches between the GOES and ER-2 lidar for the CS days. Data taken during 
October 27 was located just outside the northwestern boundary of the 1F0 
area. On October 28, the ER-2 was over Lake Michigan at - 1600 GMT, while 
it was near FMC at 1800 GMT. The three scenes at 1600 GMT correspond to 
small contiguous areas . The scattering efficiencies follow the same trends 
observed for the ground sites with the greatest values of r near local 
noon. The average value of r over the lake at 1600 GMT is 2.1, a value 
very close to that observed over land for 1600 GMT (see Table 3) . 

d. Cloud parameters from AVHRR data 

Figure 19 presents the averaged cloud emittance -albedo pairs derived from 
both AVHRR and GOES data taken over WAU on October 28. Although the 
relative viewing angles were within + 8® for the two satellites, the 
absolute viewing conditions were different. The GOES viewed the region from 
the south while the AVHRR viewed from the north. The solar plane was almost 
midway between the two orbiters. Discrepancies in the albedo range may be 
due to this misalignment (more dark pixels were found in the GOES results) . 



27 


Resolution differences , however , would tend to produce a smaller range in 
both albedo and emittance for the GOES (8 -km) relative to the AVHRR (1x4 

km^) radiances as observed in Fig. 19. In general, the GOES data are very 
close to the AVHRR results except for the smaller values of a c . The latter 

may be affected by the presence of "dark” pixels. 

Another comparison of GOES and AVHRR results is given in Fig. 20 for data 
taken over WAU November 2. The AVHRR viewed the scene with 6 - 18° and ^ 

— 38°. The squares correspond to AVHRR data as reduced for the given 
viewing angles. AVHRR emittances corrected to the GOES viewing zenith angle 
with (2) are denoted with the crosses. The range in c« c is smaller for GOES 

than for the AVHRR consistent with the resolution differences. In this 
case, the AVHRR data produce a much lower minimum cloud albedo. The 
application of (2) appears to have produced very similar emittances for the 
two data sets . 

All of the coincident AVHRR and GOES data are summarized in Table 5. The 
AVHRR Infrared optical depths are all slightly greater than the 
corresponding GOES values, while the differences in vary from scene to 

scene. Even when the times and angles are very close and the data appear 
similar as in Fig. 19, there are substantial differences in r v and r. 

Over FMC during the 28th, there is good agreement between the parameters, 
however. The outstanding differences may be attributable to a number of 
factors discussed below. 
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e. Uncertainties 

The parameter values derived here are subject to considerable uncertainty 
as evidenced by the results in Table 5 and the large standard deviations in 
earlier figures. Potential sources of error abound in an analysis of this 
type due to the large number of variables and the nonuniformity of cirrus 
clouds . 

Parameters derived from the lidar essentially provide a two-dimensional 
view of the cirrus clouds. The assumption that z c , h, and z^ represent 

the average cloud heights within the large areas covered by the strip of 

pixels is difficult to evaluate. One means of estimating how well the lidar 

data represent the large area cloud characteristics is to examine the 

differences between the strip of pixels and surrounding areas. The rms 

difference between the emittances for the strip and the box containing the 

strip is 0.05 or 7%. This difference is equivalent to a ± 0.7 km variation 

in cloud center height between the strip and the box. Changes of 2 km in 

cloud center altitude during a given half hour are common as seen in Fig. 1. 

It appears that the variations in the small scale lidar data are greater 

than those in the large scale satellite data. 

Assuming that the large scale differences are representative of the 

lidar- satellite scale differences, it is estimated that the use of lidar 

data to set z causes an uncertainty in « of + 10% based on an average 
C a 

value for e of 0.62. The accuracy will probably be less for lower, 
a. 

thicker clouds and slightly greater for high, thin clouds. The latter are 
better defined in the lidar data and present more contrast to the satellite 
view than the former. Note that no clouds with € ^0.1 were included in 
the analysis. A conservative estimate of the uncertainty in T g is + 2K. 
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Inclusion of this error raises the overall uncertainty in e & to ± 13%. 

This uncertainty in e is equivalent to a ± 20% uncertainty in r g for a 

fi 

given scene. The average IR optical thickness is 0.96 for all 71 scenes. 
From the strip and box comparison, it is also estimated that z t and h 

have uncertainties of ±0.7 km. 

Another source of uncertainty in is the use of a mean cloud height 

for the entire scene. This error source may be examined by performing a 
pixel-by-pixel analysis on a scene which varies systematically with time. 

One example is the cloud over EMC between 2020 and 2030 GMT. The GOES 
pixels from the corresponding wind strip data were averaged in lines 
perpendicular to the wind vector. Using the wind speed, these averaged 
pixels were converted to times and aligned with the lidar- defined cloud 
parameters. The results shown in Fig. 21 indicate good alignment between 
the two data sets. It appears that, in this case, the lidar data provide an 
accurate cross section of the cloud. The GOES reflectance increases as the 
cloud thickens and T increases as z c lowers. Equations (1) and (5) were 

applied to each average pixel using T £ derived from Fig. 21 to determine 

r and r . Figure 22 shows the variation of the parameters with time, 
e v 

Although the thin part of the cloud is detected with the IR data, a value 
for r is not computed since the reflectances are lower than that for 

clear skies. Nevertheless, the mean values for r & and « b derived on a 

pixel -by -pixel basis are 0.59 and 0.54, respectively, compared to 0.62 and 
0.59 derived for the entire scene using a mean value of z £ . This 

comparison suggests that the error in for using the mean cloud height 
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is around 10%. While these results may not represent all cases, they 
indicate that the use of a mean cloud height for the analysis is a 
reasonable approach . 

Assessment of the accuracy of r v relies on the satellite data and the 

models used to solve (5). The two greatest sources for error in (5) are the 
bidirectional reflectance model and the theoretical cloud albedo model. The 
latter is based on a set of calculations for one typical set of micro- 
physical properties. While the scattering phase function should not vary 
significantly for hexagonal ice crystals of various sizes, other shapes and 
phases such as bullet rosettes or water droplets may occur in the observed 
clouds. Thus, some uncertainty arises from the use of the specific model. 
Another source of uncertainty is the interpolation scheme used to extract 

T v from a set of four cu rves defining the relationship of r and a as 

v c 

a function of 8 0 . It is estimated that the use of the log-linear 
interpolation for a given 8 0 is accurate to about + 5% for the solar 
zenith angles used here. 

The anisotropy of the reflected radiation field for real clouds depends 
on the optical thickness, incident radiation, microphysical properties of 
the cloud, and the morphology of the cloud field. The value of * c used 

here is fixed for a given set of angles and represents an empirical average 
for all cloud types. Since cirrus clouds are the only type considered 
here and the angles are fixed for a given hour, it is likely that x will 

be biased with respect to local time. There will also be random errors in 
*c due to variations in microphysics, morphology, and cloud optical depth 

for a given hour. The magnitudes of these errors are currently unknown, but 
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are potentially large. If it is assumed that the time sampling represents a 
random sampling of x c> however, then averages of various parameters derived 

from all times should minimize the effects of biases in * c - 

The bidirectional reflectance model errors are estimated by reanalyzing 
the data using x Q “ 1 • The resulting overall mean value of r is 1.9, 

less than the value of 2.0 expected for large particles. The rms change in 

y of + 10% results in an 8% decrease in r and an rms difference of 18% . 
*c ~ 

Values of r near midday are increased relative to the nominal results , 
while those in the early morning and late afternoon are decreased. Since 
range in r is reduced by using the models, it may be concluded that 
their application was appropriate for this data set. 

The dispersion of the nominal ratios, however, is significant on an 
hourly basis even when the models are used. If it is assumed that r is 
Independent of time, then mean values of x c ma y t> e determined from (5) if 

r and r are known. The data were reanalyzed forcing the scattering 
e 

ratio to a value of 2.08 and solving (5) for x q - Resulting mean hourly 
values for x c are given in Table 6 with the nominal values from the 

bidirectional reflectance model. The results indicate more anisotropy in 
the cirrus reflectance pattern than in the empirical model for all cloud 
types. From these results, it may be concluded that the bidirectional 
models are accurate to within + 20%. 
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5. Discussion 

a. Infrared parameters 

The CS results in Fig. 14 are very similar to those derived by Platt 

et al. (1987; hereafter, PSD) from ground-based observations in Australia. 

The mean value of « is 0.33 for the CS , the same as that found by PSD for 

a 

their midlatitude site. The variation of vertical beam emittance with T c 

for the IFO data probably differs from that for the CS and the PSD data 
because of substantial differences in cloud depths. Clouds observed during 
the CS were 2 . 6 km thick on average compared to a mean thickness of 4.2 km 
during the remaining cirrus days of the IFO. This difference in thickness 
is evident in Fig. 15 where the CS data, a subset of the IFO data, are also 
consistent with the results of PSD. The variation of the CS volume 
absorption coefficients in Fig. 16 are not as well behaved as emittance and 
cloud thickness when considered as a function of T c< Except for the 

highest and lowest temperatures, though, for the CS data is close to 

that observed by PSD. Differences between the IFO and CS data are also seen 
in the behavior of a suggesting that thickness is not the only 

discrepancy between the CS and other IFO clouds. Cloud ice water content or 
phase may also also affect these differences. 

The well -correlated variation of r^ suggests the possibility that 

as well as T may be retrieved with VIS-IR radiance pairs. This parameter 
c 

integrates many of the other parameters examined earlier. Since cloud 

thickness is small for the highest clouds, there is little difference in 

T and T . At lower altitudes , there is more cloud thickness so it is 
c t 

possible to sense radiation from areas deep within the cloud. Whether the 
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relationship shown in Fig. 17 is typical for all observing angles is 
unknown. From this relationship and those shown in Figs. 15 and 16, it 
should be possible to obtain reasonable estimates of h and assuming 


that T is determined accurately, 
c 

The GOES IR optical depths are consistently less than the AVHRR values by 
— 0.1. This difference indicates the possibility of a calibration offset in 
the thermal channels. Despite this obvious bias, the good relative 
agreement in Fig. 20 between the corrected AVHRR emittances and the GOES 
emittances suggests that (3) is a reasonable approximation to the IR 
absorption optical depth. Any IR scattering effects which are ignored here 
are apparently insignificant compared to the other error sources. 

b. VIS parameters 

Dcrivfltion of T from a given VIS-IR radiance pair depends critically 

on the values of r and r. These parameters are subject to the model 

used in the analysis, the sizes and shapes of the particles in the cloud, 

and the cloud field geometry. The physical model used for extracting the 

cloud albedo from the measured reflectance also influences the values of 

r and r The fact that there is good agreement between the average 
v 

derived values of r over land and water at the same hour (i.e., 1600 GMT) 
suggests that (5) is a reasonable model for accounting for surface 
reflectance. Thus, other factors are most likely to cause the observed 
variations in the scattering efficiency. 

Theoretically, r approaches a value of 2.0 asymptotically from a maximum 
of — 4 when the size parameter, 2wa/A » 1 (e.g. , Hansen and Travis, 1974) . 


34 


The scattering wavelength is X , and a is the particle radius. This 
theory requires that the IR absorption efficiency is unity and no absorption 
occurs in the VIS. PSD have presented theoretical results indicating that 
r is also a function of particle shape. Though there Is considerable 
uncertainty associated with the parameters derived for a given scene, the 
averaged data presented here for both the CS and IFO scenes are consistent 
with the limit of r - 2.0. 

The results shown in Fig. 13 also indicate that r depends on T . For 

c ’ 

T c < 235 K, the average scattering ratio is r « 2.46 + 0.20, while for 

^ 235 K, r — 1.85 + 0.17. Some of this difference may be due to 

temporal sampling differences. The sampling effects may be estimated using 

the hourly mean values of r in Table 3 to compute r instead of the 

observed values. Using that approach, r - 1.96 for T > 235 K and r - 

c 

2.09 for T c < 235 K. The former value is not statistically different from 

the observed mean scattering ratio for T £ > 235 K. The latter value of r, 

however, is statistically different from 2.46, indicating that the values of 
r > 2.0 for the lower temperatures are significant. A reanalysis of the 
data using X c (r-2.08) from Table 6 also yielded temporally unbiased values 

of r for a given T c> The results using * c (jr-2.08) re tain the trend of 

increasing r with decreasing T £ (Fig. 23). This reanalysis reveals an 
almost linear trend of r v with h (Fig. 24) lending further support to the 
potential for deriving cirrus cloud thickness from the VIS-IR data. 
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Platt and Dilley (1981) found that the lidar- measured backscatter- to- 
extlnction ratio is relatively constant for T c > 233 K. The value of that 

ratio changes abruptly around 233 K to a new value. Heymsfield and Platt 
(1984) were able to show that the number density of large particles in 
cirrus clouds decreased substantially when the temperature dropped below 
-40°C. Their data also show that ice crystal habits change significantly 
when the temperature drops below 233 K. A combination of smaller crystals 
and different shapes may explain the increase in r for T c < 235K observed 

in the current data set . 

The temporal biases in r are primarily due to shortcomings in the 
analysis' treatment of X Q • Optical depth and particle shape are two of the 

more important variables determining the bidirectional reflectance. The 
average cloud optical depth in the bidirectional reflectance model used here 
is probably close to 10, while r for the clouds analyzed in this study is 

generally smaller than 2. Takano and Liou (1989b) indicate that the 
reflectance anisotropy diminishes with increasing optical depth for 
scattering by hexagonal columns . This effect may help explain the larger 
range of x c ( r “2.08) compared to the nominal values of x c in Table 6. 

The nearly monotonic increase of x c with 8 in Table 6 for both cases 

is consistent with the scattering phase function for randomly oriented 
hexagonal columns reported by Takano and Liou (1989a) . Since multiple 
scattering effects tend to smooth out some of the smaller variations found 
in the single -scattering phase function, a direct correspondence between the 
phase function and x c is not expected. The prominent trends, however, 

should be apparent for the optical depths considered here. The minimum 
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value of x at 8 - 120° in the afternoon is near the minimum in the phase 

function of Takano and Liou (1989a), while the distinct maximum near noon 

corresponds to the backscatter relative maximum in the theoretical model. 

Mie calculations for a sphere with a 10-/im radius would produce an increase 

in y with 8 between 100° and - 135° followed by a flattening of the 
c 

values of x c for greater scattering angles until 8 > 170° where the glory 

becomes apparent (Hansen and Travis, 1974). These trends suggest that a 
model which simulates reflected radiances for a cirrus cloud using a phase 
function for spheres may bias the values of derived from measurements 

taken at a given 8. 
c. "Dark" pixels 

It is possible that the "dark" pixels noted earlier may be caused by the 
extreme reflectance anisotropy expected for very thin clouds. For small 
values of r and 8 - 120°, there should be very little reflection from 

the cloud in the direction of the satellite compared to its albedo to 
upwelling diffuse radiation from the surface. Thus, the pixel may appear 
darker than a clear scene for the same viewing and illumination conditions . 

Although such a mechanism may cause some pixel darkening, it is probably 
not the primary one. Cloud shading is more likely to be the predominant 
cause of dark pixels. For example, at altitudes of 10 km, a 2-km thick 
cloud can cast a shadow which has its leading edge removed at least 8 -km 
horizontally on the surface from the position of the cloud's leading edge 
for 6 0 > 45°. Thus, the cloud and the portion of the surface affected by 

direct transmission through the cloud can easily be offset by one or more 
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GOES pixels. The projected line of site from the satellite through a cloud 
will not line up with the cloud's shadow unless 6 » 180°. In most 
instances, therefore, the direct transmission from the surface through the 
viewed cloud will come from the shadow of a different cloud. The diffuse 
reflectance will result from the combined effects of the various clouds in 
the vicinity. For a homogeneous cloud field, an optically thick viewed 
cloud, or 6 - 180°, such effects are negligible. Figures 1, 2, and 21 
provide ample evidence that the clouds observed during the CS are neither 
optically thick nor homogeneous. In an inhomogeneous cloud field, it is 
possible to observe a surface darkened by the shadow of a thick cloud 
through a thin cloud which produces little scattering of its own in the 
direction of the satellite. Thus, to the viewer in space a pixel appears to 
be a cloud by virtue of its cold temperature, but is darker than expected 
for a clear scene. 

This effect may be examined quantitatively by considering (5). Use of 
that model implicitly assumes a homogeneous cloud field. To consider an 
inhomogeneous cloud field, let t replace in (6) , where r s is the 

VIS optical depth of the cloud casting a shadow onto the surface in the line 
of site of the observed cloud which has optical depth r . Since the 

diffuse radiation emerging from the bottom of the clouds and reflecting from 
the surface is coming from all of the clouds in the field, let a c in the 

third term of (5) be replaced by the average albedo of the observed and the 
shading cloud. The results of solving this inhomogeneous version of (5) for 
two different observed clouds are displayed in Fig. 25. The circles 
represent r g - r^. Viewing angles were fixed at 6 - 52° and x c varied 
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with 9 0 as in the lower half of Table 3. Clear-sky reflectance is denoted 
with the dashed horizontal line. 

The greatest effects of cloud inhomogeneities are seen for - 0.25. 

Both shading clouds cause dark pixels for 6 0 < 81°. Some dark pixels occur 

for r - r ^ for 9 0 < 60°. Calculations using smaller p 0 ' s produced no 

dark pixels. The thickness of the shading cloud becomes more important as 
fi 0 increases . The impact of shading on the observed cloud having - 1 

is less pronounced with little likelihood of dark pixels for the range of 
H 0 considered. Shading by the thin cloud (r - 0.25) actually increases 

the observed reflectance for all cases using r - 1. As increases to 

larger values, the impact of r g will become negligible. This lack of 

shadowing effect for larger clouds suggests that the observed reflectance 
may be biased toward a value which is lower than expected when the scene 
contains both optically thick (r > 2) and thin clouds, even if no dark 

pixels occur. 

Shadowing effects can be observed in the data shown in Fig. 22. Dark 
pixels corresponding to the thinnest part of the cloud are seen before 2027 
GMT. The next few pixels to the right are detectable but r is much less 

than r . The VIS optical depth for the pixels corresponding 2040 GMT is 

more than three times the value of r suggesting enhanced reflectance. In 

this figure, the sun would be located approximately to the right of the 
figure at 8 0 - 68°. With a cloud top near 11 km, the thickest part of the 

cloud would cast a shadow approximately 28 km to the left or to the position 
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corresponding to 2025 GMT. The view from the satellite is at a 35° angle, 
not perpendicular, to the solar plane. Thus, the surface in the line of 
site of the cloud at 2025 GMT would correspond to the 2020 GMT position. 

That surface would be shaded by a cloud with r g « 0.5 or r v « 1.0 assuming 

that Fig. 22 is an accurate cross-section of the cloud. A dark pixel would 

be expected at 2020 GMT based on the results in Fig. 25. The first visibly 

detectable cloud coincides with r - 0.3 at 2027 GMT. Its value of r 

e v 

is much less than the expected value of 0.6. 

While the shading can explain the darker pixels, it does not account for 

the excessively bright ones at 2040 GMT. For these bright pixels, it is 

instructive to refer to the cloud structure in Fig. 1. The vertical 

thickness of the cloud changes rapidly so that the pathlength of the 

incident solar radiation through the cloud is not necessarily hsec0 o . The 

pathlength may be substantially increased because of the cloud structure. 

For example, the sun angle may coincide with the right arm of the "V" 
defining the cloud centered at 2035 GMT in Fig. 1. This effect would tend 
to cause a cloud to appear brighter than expected for the observed r . In 

this case, the shadowing and enhancement effects apparently cancel to yield 
a reasonable value of r for the whole scene. 

The occurrence of dark pixels is shown in Fig. 26 as a function of 

Most of the dark pixels are found at emittances corresponding to < 0.2. 

Some, however, are found for r - 0.5. Very few dark pixels were observed 

near noon when 9 > 160°. Most were found in the midmorning and late 

i 

afternoon when shading conditions were favorable. The dark pixels comprise 
only 3% of the data considered here. This percentage belies the importance 
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of this effect since shading will occur in many instances without producing 
the easily detectable dark pixels. The shading effect will tend to reduce 
the observed reflectance causing an underestimation of the cloud s thickness 
and emittance . Cloud shadows or their absence will also affect the 
interpretation of reflectance anisotropy. Their presence at high solar 
zenith angles will yield lower than normal values of r y . 

6. Concluding Renarks 

The cirrus clouds observed during the case study days are similar to 
those observed in previous research. By comparison, however, they are not 
necessarily typical of midlatitude cirrus clouds since their bulk properties 
are significantly different from those observed during other IFO cirrus 
days . 

The analysis performed in this study using a combination of ground-based 
and satellite instrumentation yielded some results which are similar to 
earlier studies using other combinations of platforms and instruments. This 
consistency of results for different approaches lends a higher degree of 
confidence to the common findings. In addition, the combined data sets used 
here have provided some valuable new insight into the problem of determining 
cirrus cloud properties using VIS and IR radiance pairs from satellites. 

The determination of cloud- top emittance or actual cloud- top temperature 
appears to be feasible based on the results of this study. Previously, what 
has been termed cloud-top temperature in emittance -adjusted VIS-IR 
retrievals is closer to to cloud center temperature. The emittance ratio 
and the dependency of cloud thickness on cloud-center temperature derived 
here may be utilized in a scheme to estimate cloud thickness. This latter 
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parameter will be useful for computing radiative divergence for observed 
cirrus clouds . 

The results indicate that scattering efficiencies are greater for colder 
clouds. This finding holds promise for improvements in determining the IR 
optical depth from reflected VIS data. Additional study, however, is 
required to confirm this conclusion. 

From the examination of the reflectance data and scattering efficiencies 
it is concluded that much work remains to adequately describe the scattered 
radiation field for real cirrus clouds. Both theoretical and empirical 
bidirectional reflectance models should be developed for clouds composed of 
realistic particle shapes. As in the ISCCP model, the patterns should be 
developed for various optical depths. Those models will require 
confirmation with the aid of further observational data. The effects of 
cloud particle scattering in IR radiative transfer should also be examined 
observationally . 

Cloud shadows are a problem for the interpretation of cloud reflectances, 
especially for cirrus. The effects are not limited to large solar zenith 
angles. Relatively high viewing zenith angles can produce situations which 
cause the viewing of shadows even for the near-zenith sun. It is apparent 
that analysis of a single pixel is most subject to shadowing problems. Some 
of the effects may be diminished through averaging over several pixels. 

Other approaches to dealing with shading need to be developed. Validations 
of cirrus scattering calculations must also consider these effects. The 
problem of shadows is not as important over water because of the low surface 
albedo. Other problems in remote sensing of clouds such as partially cloud- 
filled pixels have not been considered here. Future research efforts should 
address these other factors. 
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APPENDIX 


List of Symbols 


a 

B 


h 

k, k 0 


k 


max 

r, r 

T p T 


e 


s 



T 

max 

t 



£ b 


particle radius 

Planck function at 11.5 pm 

observed and clear -sky visible counts 

cloud thickness 

albedo coefficient, ozone absorption coefficient 
albedo coefficient for maximum high-cloud albedo 
scattering efficiency ratio, cloud- top/-center emittance ratio 
observed and clear-sky 11.5 pm equivalent blackbody temperatures 
cloud-center and -top 11.5 pm equivalent blackbody temperatures 
clear-sky temperatures estimated from surface and satellite data 
air temperatures at altitude z and at the surface 
maximum observed 11.5 pm equivalent blackbody temperature 
time 

ozone abundance 

cloud-base, cloud-center, and cloud- top altitudes 
cloud and diffuse cloud albedos 
clear-sky albedo, clear-sky diffuse albedo 
maximum albedo for high clouds at a given emittance 
volume backscatter coefficient 

vertical effective emittance, effective beam emittance 


6 


single -scattering angle 
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e, e 0 


viewing zenith angle , solar zenith angle 


A latitude or wavelength 

p, p 0 cos 6, cos8 0 


p ' p s' p c 


o 

a 


r 


e * 


r 


v 


T 

S 


observed, clear- sky, and cloud visible reflectance 
volume extinction coefficient 

infrared absorption and visible scattering optical depths 
visible scattering optical depth of shading cloud 



ozone, downward cloud, and upward cloud transmittances 


* 



CS 

IFO 

IR 

FMC 

GMT 

MAD 

WAU 

VIS 


longitude 

clear- sky and cloud anisotropic reflectance factors 

relative azimuth angle 

Case Study (October 27-28, 1986) 

Intensive Field Observations (October 19 - November 2, 1986) 
infrared (11.5 /im) 

Ft. McCoy, Wisconsin 
Greenwich Mean Time 
Madison, Wisconsin 
Wausau, Wisconsin 
visible (0.65 pirn) 
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Table 1. Times and locations of lidar-satellite data used in this study. 


Site 

Month 

Day 

Times (GMT) 

FMC 

October 

22 

1300, 1330, 1400, 1600, 1630, 1700, 2000 



27 

2030, 2100 



28 

1330, 1400, 1430, 1500, 1600, 1700, 1900 




1930, 2000, 2030, 2100, 2130, 2200 



30 

2000, 2030 


November 

1 

1800, 1900 



2 

1900, 2000, 2100 

MAD 

October 

28 

1330, 1500, 1600, 1700, 1800, 1930, 2000 




2030, 2100, 2130, 2200 

WAU 

October 

22 

1300, 1330, 1400, 1430, 1600, 1630, 1700, 1800 




1830, 1900, 1930, 2130, 2200 



28 

1500, 1600, 1700, 1800, 1900, 1930, 2000, 2030 




2100, 2130 



30 

2130, 2200 


November 

1 

1800, 1900, 2000 



2 

1700, 1800, 1900, 2000, 2100 











Table 2. Observed and computed cloud properties for October 1986 Case Study. 
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1600 

283.0 

10.0 
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1930 

285.5 
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276.8 
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1800 
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1930 
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Table 3. Reflectance parameters computed for all data. 
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1400 

1430 

1500 

1600 
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0.83 

0.28 

0.68 
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Table 4. Observed and computed cloud properties for Case Study ER-2 data. 



Table 5. Comparison of AVHRR- and GOES -derived cloud parameters. 
































Table 6. Anisotropic reflectance factor comparison. 
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Outline of cirrus cloud derived from lidar baclcscatter intensities 
over FMC , October 28, 1986. 

Estimated lidar volume extinction efficiency from cirrus cloud 
over WAU, October 28, 1986. 

Comparison of clear- sky and shelter air temperatures over FMC, 
October 27, 1986. 

Theoretical observed VIS reflectance over WAU for 1A00 GMT 
(p - 0.13) and for 1800 GMT (p o - 0.16). 

Theoretical observed albedos for various surface albedos. 

Example of low-cloud filter application over FMC at 1930 GMT, 
October 28, 1986. 

Clear histogram over FMC at 1500 GMT, October 27, 1986 (numbers 
denote frequency of occurrence of temperature -count pairs) . 
Histogram of cirrus clouds over FMC at 1500 GMT, October 28, 1986. 
Cloud emittances and observed reflectances for T < T g • 3 K 
derived from Fig. 6b. 

Cloud emittances and cloud albedos derived from Fig. 6b without 
"dark” pixels. 

Average cloud emittance as a function of cloud albedo derived 
from Fig. 7b. 

Average cloud emittance versus cloud albedo at 1500 GMT, October 
28, 1986 over three sites. 

Average cloud emittance and albedo at three times during October 
28, 1986 over FMC. 

Average cloud emittances and albedos at 2100 GMT over FMC during 
different days. 
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Fig. lib. 
Fig. 12. 
Fig. 13. 
Fig. 14. 

Fig. 15. 
Fig. 16. 

Fig. 17. 
Fig. 18. 

Fig. 19. 

Fig. 20. 
Fig. 21. 

Fig. 22. 
Fig. 23. 
Fig. 24. 

Fig. 25. 

Fig. 26. 


Average cloud emittance versus cloud albedo for all case study 
scenes, T z - T c< 

Same as Fig. lib, except for all IFO scenes. 

Variation of mean scattering efficiency with IR optical depth. 
Mean scattering efficiency versus cloud-center temperature. 
Variation of cloud vertical beam emittance with cloud- center 
temperature. 

Dependence of cloud thickness on cloud-center temperature. 

Cirrus volume absorption coefficient versus cloud-center 
temperature . 

Same as Fig. lib, except T z - T c> 

Dependence of ratio of cloud- top emittance to cloud- center 
emittance on cloud-center temperature. 

Comparison of cloud albedos and emittances derived from October 
28, 1986 GOES and AVHRR data taken at - 2100 GMT over WAU. 

Same as Fig. 19, except for November 2, 1986 at - 2000 GMT. 
Comparison of GOES and lidar observations along wind vector over 
FMC during October 28, 1986. 

Cloud optical properties derived from Fig. 21. 

Same as Fig. 13, except for x c (r-2.08). 

Variation of cloud thickness with VIS optical depth for 
* c (r-2.08). 

Effect of cloud shading on observed reflectance. Solid symbols, 
r — 1. Open symbols, r v - 0.25. 

Percentage of total number of "dark" pixels as a function of 
effective beam emittance, 
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